arXiv:1506.00550vl [astro-ph.CO] 1 Jun 2015 


A novel approach to reconstructing signals of isotropy violation 

from a masked CMB sky 

Pavan K. Aluri,^ Nidhi Pant/ Aditya Rotti/’^ and Tarun Souradeep^ 

^lUCAA, Pune - 411007, India 
^Florida State University, Tallahassee, FL 32304, USA 

Statistical isotropy (SI) is one of the fundamental assumptions made in cosmological model build¬ 
ing. This assumption is now being rigorously tested using the almost full sky measurements of the 
CMB anisotropies. A major hurdle in any such analysis is to handle the large biases induced due to 
the process of masking. We have developed a new method of analysis, using the bipolar spherical 
harmonic basis functions, in which we semi-analytically evaluate the modifications to SI violation 
induced by the mask. The method developed here is generic and can be potentially used to search 
for any arbitrary form of SI violation. We specifically demonstrate the working of this method by 
recovering the Doppler boost signal from a set of simulated, masked CMB skies. 
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I. INTRODUCTION 

Cosmic microwave background (CMB) observations 
have ushered in the precision era in cosmology. This 
enables us to test the basic premises on which the stan¬ 
dard cosmology rests, such as the assumption of statisti¬ 
cal isotropy via the Cosmological Principle. Some claims 
of statistical isotropy (SI) breakdown, that appeared in 
CMB literature were revisited by WMAP and PLANCK 
satellite science teams; see Ref. [1, 2] and the references 
there in. 

The Bipolar Spherical Harmonic (BipoSH) formalism 
[3, 4] is a powerful tool to identify and study statistical 
isotropy (SI) violation. This has been one of the methods 
used to search for deviation from isotropy in the CMB 
data [1, 2, 5]. 

The idea behind searching for SI violation is simple : it 
involves searching for non-vanishing BipoSH coefficients. 
This is made non-trivial due to the presence of a mask, 
as it induces spurious correlations that aliases the mea¬ 
surements. In this article, we present a new methodology 
based on the BipoSH formalism, to estimate signals of SI 
violation from a masked CMB sky. The algorithm pre¬ 
sented here is generic, and can be potentially applied to 
recover any signal of SI violation, such as weak lensing, 
hemispherical power asymmetry etc. The Doppler boost 
is a known isotropy violation in the CMB at a challenging 
level of subtlety. In this article we demonstrate the re¬ 
covery of the Doppler boost field from a Doppler boosted 
CMB sky in the presence of a mask using our formalism. 

The paper is structured as follows. We briefly recap 
the BipoSH formalism in Sec. H. In Sec. HI we first dis¬ 
cuss the SI violation induced due to Doppler boosting 
and masking separately and finally the resultant SI vio¬ 
lation when both the effects are considered together. In 
section IV we discuss the minimum variance estimator 
that takes into account the effects of masking for an ar¬ 
bitrary signal. The performance of the new estimator is 
discussed in Sec. V, with the specific example of Doppler 
boost. Finally our conclusions are presented in Sec. VI. 


II. BIPOSH : A BRIEF SUMMARY 

A CMB anisotropy map, T{n), defined on a sphere is 
conventionally decomposed in terms of spherical harmon¬ 
ics (Yi^) as, 

OO +/ 

T{n) = '^ X , (1) 

1 — 1 m——l 

where h denotes the position on the sphere, T{h) are the 
temperature anisotropies observed in the direction fi and 
aim are the spherical harmonic coefficients of expansion. 
The two-point correlation function is defined as, 

C{hi,n2) = {T{ni)T{n2)), ( 2 ) 

where (■ ■ •) denotes an in principle average over an en¬ 
semble of statistically independent CMB realizations. If 
the field T is statistically isotropic, then it can be argued 
that the correlation function cannot have any explicit di¬ 
rectional dependence, and that it can only be a function 
of the separation 9 between the two directions hi & h 2 , 
where cos 0 = hi • h 2 . As a result of this simplification, 
the two-point correlation function is given by, 

I OO 

C(hi,h2) = C(0) = — ^(2Z + l)GPi(cos0), (3) 
1=1 

where Pi are the Legendre polynomials and Ci is the well 
known angular power spectrum. When expressed in terms 
of the harmonic space covariance matrix, 

~ OJlSii' 6mm' , (4) 

Cl ’s are the diagonal and the only non-vanishing elements 
of the matrix. Hence for a statistically isotropic Gaussian 
random field, C/ completely characterize the statistical 
properties of the field. 

In the absence of statistical isotropy, Ci do not com¬ 
pletely describe the CMB sky and the Bipolar spherical 
harmonic (BipoSH) basis functions [3] are better suited 
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for this case. BipoSH form a complete orthonormal basis 
for an x 5^ space, and hence can be used to describe 
any two point correlation function for a field defined on 
a sphere as 

C{ni,n2)= ^hh{^hini)<^Yi2in2)}LM, ( 5 ) 

L,M,h,h 

without making assumptions about statistical isotropy 
of the field. Here are the coefficients of expansion 
and {Yi,(^YiJlm = are the 

bipolar spherical harmonic basis functions [6]. 
are the Clebsch-Gordon coefficients whose indices satisfy 
the properties : a) \li — h] < L < \li +Z 2 I, b) toi +TO 2 = 
M, c) —li < mi < li, and d) —I 2 < m 2 < h- 

The harmonic space covariance matrix for a statisti¬ 
cally anisotropic field, is fully described in terms of the 
BipoSH coefficients as, 

{0-1117120-121^12) = ^ iimii 2 m 2 J (®) 

LM 

in complete analogy with Eq. [4]. Therefore the BipoSH 
coefficients completely characterize the statistical proper¬ 
ties of a Gaussian, statistically anisotropic random field. 
It can be shown that = (—1)*-\/2/ -I- IG/. Hence the 
BipoSH coefficients are a generalization of the commonly 
studied angular power spectrum. Note that Eq. [6] can 
also be used to define an estimator for the BipoSH coef¬ 
ficients as follows, 

aLA4 \ '' /^LIvI ('y\ 

^hh ~ / > ^Ii7nil2m2^h'i-mi0il2m2 ' 

mi?7i2 

It is easy to see that this is an unbiased estimator, since 
it converges to the true BipoSH coefficients on averaging 
over an ensemble of CMB skies. These BipoSH coeffi¬ 
cients have the following symmetry properties, 

• Exchange symmetry : G(ni,n 2 ) = G(n 2 ,ni) 

( 8 ) 


• Reality condition : C*(hi, 712 ) = C(hi,h 2 ) 


A 


LM 
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(_l)Af(_i)/i-Hi2-L^ 


L,-M 

I 1 I 2 


(9) 


Alternatively, following the triangularity condition be¬ 
tween {li,l 2 ,L}, it is convenient to define the BipoSH 
coefficients as where we set h = I, and I 2 = l + D, 

with the constraint that \D\ < L. Needless to say, fol¬ 
lowing the above mentioned properties of the BipoSH 
coefficients, determining the BipoSH coefficients for 0 < 
D < L and 0 < M < L will amount to computing all the 
BipoSH coefficients, 

BipoSH have been a powerful tool in exploring isotropy 
violation in a blind as well as honed manner [7]. It has 
earlier been used to probe SI violation as induced due to 
specific models of isotropy violation (phenomenological), 
cosmic topology (cosmological) and beams (systematic) 
[ 1 , 2 , 8 - 10 ]. 


III. EFFECT OF ISOTROPY VIOLATING 
PHENOMENA 

Here we analytically derive the form of BipoSH coef¬ 
ficients as a function of multipole I, induced by Doppler 
boosting of the GMB sky. Following this, we derive the 
effect of masking on BipoSH coefficients, in addition to 
the anisotropy due to Doppler boost. Based on these an¬ 
alytic studies, a new estimator is obtained to recover the 
Doppler boost vector from a masked GMB sky. 


A. BipoSH coefficients due to Doppler boosting of 
CMB anisotropies 


One of the known phenomena that leads to break¬ 
down of statistical isotropy is our relative motion (v) 
with respect to the CMB rest frame (i.e. the frame in 
which an observer does not see a dipole). The strongest 
anisotropy in an otherwise uniform background temper¬ 
ature of To = 2.7255 Kelvin [11, 12], is an excess (deficit) 
temperature in the direction (opposite direction) of this 
relative motion. Doppler boost (j3 = v/c) that induces 
this large scale dipole temperature anisotropy (=Tol3 ■ h) 
has a well measured amplitude and direction given by 
|/3| = 1.23 X 10“^ and /3 = (264°,48°) in galactic co¬ 
ordinates, respectively [13, 14]. In addition to this dipole 
anisotropy, the velocity boost also modulates and aber¬ 
rates the CMB temperature fluctuations leading to SI 
violation [15-20]. These aberration and modulation ef¬ 
fects can be used to make an independent measurement 
of the velocity of our local motion. It has been measured 
with ESA’s PLANCK satellite data also [21, 22]. 

Doppler boosted CMB temperature anisotropies are 
given by [21], 

T(h) = 7o/ 3 • n -I- 0 — V(/3 • h)^ ^1 -I- ■ h^ , (10) 


where T and 0 correspond to boosted and unboosted 
temperature anisotropies respectively, and f3 = v/c 
is the Doppler boost vector. From Eq. [10] we see 
that, in addition to the dominant dipole component due 
to the Doppler effect, Doppler boosting also aberrates 
and modulates the CMB anisotropies. The modula¬ 
tion component has a frequency dependent factor hn = 
(v/yo) coth(i2/2r'o) — 1, where vq = ksTo/h « 57 GHz. 

For the rest of the discussion we will not bother with 
the first term in Eq. [10] since we often work with dipole 
subtracted maps. Evaluating the two point correlation 
function for the temperature anisotropies of a Doppler 
boosted CMB sky described in Eq. [10] and ignoring the 
dipole term, it can be shown that the resultant BipoSH 
coefficients can be cast in the following form. 


aLM 


i^hh)ub.cmb + , 


( 11 ) 


up to first order in the isotropy violating Doppler boost 
field, /3(h) = j3 ■ h. Note that Doppler boosting gen¬ 
erates only the L = 1 BipoSH modes. The term 
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(A^¥) ^ ^ corresponds to BipoSH coefficients due to 

unboosted CMB anisotropies, Plm the spherical har- 
monic coefficients of /3(n), and is the shape function 
corresponding to Doppler boosting given by, 


^hh - {^hl 2 )mod i^hh)abr 

where, 

(^^^ 2 ) mod 
(^^ 1 ^ 2 ) abr 


Gl + Cl^ 11/^11/2 


Z 10/20 ’ 


•s/Itt IIl 
[Cij^F(li, L, I 2 ) + Ci 2 F{l 2 , L, ^i)] 


( 12 ) 


(13a) 


n, 


\/47T 

Qrobo> (13b) 


F{k,Lj 2 ) = + + + + 


where Ui = ^21 + 1 and, and are 

the shape functions due to modulation and aberration 
effects, respectively and are completely specified by the 
angular power spectrum G;. 

A full sky minimum variance estimator for the veloc¬ 
ity boost can then be defined as a weighted linear com¬ 
bination of the observed BipoSH coefficients, See 

Appendix A for details. 


B. Effect of masking on BipoSH coefficients 


the CMB sky. This equation is a generalization of the 
MASTER kernel [23] which describes how the angular 
power spectrum (G/) is modified by mask. In the special 
case of L = 0 and L' = 0, this equation exactly reduces to 
the MASTER equation relating full sky power spectrum, 
G/, to partial sky power spectrum, G;. 


C. BipoSH of a masked anisotropic sky 


The full sky BipoSH coefficients of an anisotropic 
Doppler boosted CMB sky are given by Eq. [11]. Thus 
the BipoSH coefficients of a masked Doppler boosted 
CMB sky can be obtained by substituting Eq. [11] in 
Eq. [15] to get. 


= UtfO , + E l3L'M'Kt'Xi2 ^ ( 16 ) 

\ / ub.cmb 

L'M' 

where, 

\ /ub.cmb v 47 r 

h 

\ ' H^H^e/^iiO r^hO jxrLM J ^5 ^6 L\ 

^ 2 ^ ^ ^hOhO^hOleO^hh \l 2 h I 3 j ^ 1 G 


Even though an ideal CMB sky is isotropic, masking 
the sky to avoid foregrounds makes it highly anisotropic. 
In this section we derive the BipoSH coefficients arising 
due to the process of masking. Let a masked CMB sky 
be defined as, 

f (n) = W(h)r(h), (14) 


denotes the BipoSH coefficients generated due to masking 
an unboosted CMB sky characterized by G;. /3 l'M' are 
the harmonic coefficients of the Doppler field /3(h) = /3-h, 
and the term {} 2 x 3 denotes the 6j—symbol. The modi¬ 
fied shape function (MSF), of the Doppler sig¬ 

nal, is the masked analogue of the full sky shape function, 
™ [12] and is given by. 


where f and T denote the masked and unmasked CMB 
sky respectively and W(h) denotes the mask used for the 
analysis. Evaluating the two point correlation function 
for the masked CMB sky, it can be shown that the resul¬ 
tant BipoSH coefficients are given by. 


A: 


LM 

hh 


E n^Hu \" n^Hie wio W20 

^3^4 ^5^6 





(15) 


where denotes the BipoSH coefficient for masked 

CMB sky, Af'^iY denote the BipoSH coefficients of the 
full sky map, which is not assumed to be statistically 
isotropic for generality, are the BipoSH coefficients 

of the mask and { jsxs denotes the 9j—symbol. The 
mask BipoSH coefficients are defined in the same way as 




n^iiu Ubiiie f-iho 

V4)r 


I 3 I 4 I^Iq 

( L h I 2 ] 

X^UO°eoE (18) 


JK 


J I 3 Iq 


The modified shape function K’j^ 2 ii[i 2 incorporates the 
mixing of modes {JK} due to mask and the intrinsic 
anisotropic modes {L'M'}, giving rise to the observed 
modes {LMj. This coupling is captured by the Clebsch- 
Gordan coefficient and the 9j—symbol. 

Hence from Eq. [16], it can be seen that masking may 
leak power from an intrinsic anisotropic mode, L', to any 
observed mode, L. Further note that the modified shape 
function due to masking is now phase (M,M') dependent 
unlike in Eq. [12]. Here we emphasize that the formalism 
discussed in this section is completely generic, since we 
assumed no particular form of isotropy violation. 
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IV. PARTIAL SKY ESTIMATOR 
A. Approximations 

In general, masking leads to highly entangled modes as 
seen from Eq. [16] and Eq. [18]. We now specifically con¬ 
sider the case of BipoSH coefficients generated due to a 
Doppler boosted CMB sky. Though in a full sky Doppler 
boosted CMB sky the signal is only in i' = 1 mode of 
BipoSH coefficients, the masked Doppler boosted sky can 
have power leaked to L ^ 1 modes. The mode coupling 
of intrinsic anisotropic signal {L' = 1) with the mask 
indices (J), giving rise to the observed modes (L) are 
tabulated in Table I. 


Intrinsic anisotropic signal, L' — 1 


Observed BipoSH index, 
L 

Mask BipoSH index, 
J= jL-L'j toL + L’ 

L = 1 

J = 0,1,2 

L = 2 

J = 1,2,3 

L = 3 

J = 2,3,4 


TABLE I: Illustrated here are the observed BipoSH indices 
(L) that result from the mixing of an intrinsically anisotropic 
signal in L' = 1 with the mask BipoSH indices (J). For 
a mask which is largely azimuthally oriented, the dominant 
mask modes contributing to the mixing kernel are highlighted 
in bold letters. 

We now argue, that even though there is coupling be¬ 
tween different modes, it is reasonable to assume that 
only modes with L' = L = 1 and M = M', which we 
refer to as the diagonal approximation, are sufficient to 
evaluate Eq. [16]. 

Most masks used in CMB analysis can at first order 
be approximated as a band along the galactic equator 
where the astrophysical emission from our own galaxy is 
the highest and hence have a significant azimuthal sym¬ 
metry. It can be argued that for such a mask, the mix¬ 
ing kernel is primarily diagonal in {M, M'}, motivated 
by the fact that for a band mask, which has perfect az¬ 
imuthal symmetry, the {M, M'} coupling can be shown 
to be identically zero. For mask BipoSH coefficients it is 
seen that the J — 0 mode is the most dominant mode. 
Observing the various couplings in Table I it can be seen 
that only the L = 1 mode couples to L' = 1 mode via 
J = 0, while all higher L > 1 modes couple to L' = 1 
via J > 0, suggesting that the contribution from terms 
where L ^ L' will be sub-dominant. 

We follow these heuristic arguments by quantitatively 
showing that the off-diagonal terms are infact subdomi¬ 
nant by explicitly evaluating the kernel . In Fig [1] 

we compare the full sky shape function of Doppler 
boosting with the modified shape function, , eval¬ 

uated using the apodized version of the mask shown in 
Fig. [2] (used for all analyses presented in this article). 
Note that the diagonal terms are nearly two orders of 



Multipole. I 

FIG. 1: In this figure we compare the full sky shape function 
and the modified shape function modified 

due to mask, for the Doppler boosted CMB sky. The full sky 
shape function is independent of phase as it does not depend 
on M. The most dominant terms of the modified shape func¬ 
tion is the diagonal denoted by i.e. when L = L' 

and M = M' . Also plotted are some of the off-diagonal terms 
of the shape function, i.e. L ^ L' and M yf M', which 
are seen to be negligible when compared to the diagonal and 
hence ignored in the analysis. Further note that for the di¬ 
agonal terms, there is a M dependent suppression of power 
as compared to The varying supression of power be¬ 

tween modes corresponding to different M is determined by 
the morphology of the mask. 

magnitude larger than the off-diagonal terms. This jus¬ 
tifies the approximations suggested above. 

Here it is also important to note that the modified 
shape function is different for different phase modes 
(M = 0,1). This distinction is purely determined by the 
details of the mask. As we will see in the following sec¬ 
tions, this difference in MSF corresponding to different 
phase modes (M), naturally compensates for the differ¬ 
ential loss of signal in these modes while recovering the 
Doppler field harmonics. 

Note that in Eq. [16], it is fairly complicated to al¬ 
gebraically solve for the Doppler field harmonics, /3 lm, 
since it appears inside a summation. By making the ap¬ 
proximations discussed in this section, we have simplified 
Eq. [16] to the form, 

( 19 ) 

\ / ub.cmo 

This is particularly useful, since it can be used to define 
a Doppler estimator, similar to the full sky estimator. 


B. Weighted variance estimator using modified 
shape function (WV-MSFE) 

In this section we briefly discuss the Doppler estimator. 
Eq. [19] can be inverted to arrive at an estimator for the 
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Doppler field harmonics, 


Plm = 


alm _ / alm\ 

^hh /■u.b.cmb 


t^LM 

^LMh 


( 20 ) 


where denotes the BipoSH coefficients derived from 
the data maps, and (A;^^)ub.cmb is the expected bias due 
to mask and spatially varying noise in an observed map. 
The expected bias is estimated from an ensemble of sim¬ 
ulations which are not Doppler boosted. This naive esti¬ 
mator can be optimized by minimizing its variance aris¬ 
ing due to cosmic variance and instrument noise. The 
derivation of the estimator is discussed in Appendix B. 
The estimator used to reconstruct the Doppler boost vec¬ 
tor is given by. 


^LM 

PlM = wY , ( 21 ) 

where ub.cmb are the bias corrected 

BipoSH coefficients. The weights which minimize 
the variance are given by the expression. 


Kh = 


E 


^LM 


M \ 

' ^LMlil2 


-\ -1 


E- 


(_dA 


M I , , 

V LM'lil' 


where. 


! _ /I aLM \2\ I / ^LM\ 

K'^hh / “ I /ub.cmb \\-^l2l2 )ub.cmb 


( 22 ) 

(23) 


is the variance of unboosted map’s BipoSH coefficients. 

We see that the effective weights are M dependent ow¬ 
ing to the M dependent shape function in Eq. [21], unlike 
in the full sky estimator. Since the estimator is weighed 
by the modified shape function which accounts for loss of 
power due to mask, it naturally corrects for the reduced 
amplitude of the Doppler vector due to masking. 


V. DEMONSTRATION OF THE WV-MSF 
ESTIMATOR 


We demonstrate the working of this newly proposed 
method for reconstructing the Doppler boost, by evalu¬ 
ating it on an ensemble of Doppler boosted simulations. 
We also run a parallel analysis on full skies to allow direct 
comparison of it’s efficiency. 


A. Generating simulations 

We generated a set of 1000 simulations with character¬ 
istics of Planck 217 GHz instrument using the best fit the¬ 
oretical Cl corresponding to the cosmological parameters 


from Planck 2013 data [24]. The isotropic simulations 
were generated using the synfast facility of HEALPix^ 
[25]. The Doppler boosted simulations are generated us¬ 
ing the Code for Non-Isotropic Gaussian Sky (CoNIGS) 
algorithm [26]. We injected Doppler boost with ampli¬ 
tude |/1| = 1.23 X 10“^, pointing towards the galactic 
coordinates {l,b) = (264°,48°). Since these simulations 
are specific to the 217 GHz channel, we set the frequency 
dependent modulation factor to = 3. 

The noise simulations were generated separately and 
added to the GMB simulations. In our analysis, we 
first worked with isotropic noise approximation and then 
with realistic spatially varying noise, so as to progres¬ 
sively increase the complexity of the data. The isotropic 
noise simulations were generated using a white Gaussian 
noise spectrum corresponding to the 217 GHz instrument 
(Of whm = 5 arcmin & cr = 4.8 fiK/K), as quoted in 
Planck Bluebook [27]. To simulate the full mission (five 
surveys) noise levels, we divide the nominal (two surveys) 
noise standard deviation by y/S/S. 

The realistic, spatially varying noise maps were gener¬ 
ated using the noise variance map (no cross-pixel correla¬ 
tions) corresponding to Planck 217 GHz channel available 
in the public archives^. 


B. Mask and modified shape function 

We use the common analysis mask used in the Planck 
analysis shown in Fig. [2]. Sharp {0,1} masks are not 
ideal for any band limited analysis as they result in heavy 
ringing at the edges of the mask. To avoid this unneces¬ 
sary complication we apodize the mask with a Gaussian 
beam of 30 arcmin. The apodized mask has an effective 
sky fraction of fsky ~ 0.78 which is practically same as 
the sky fraction available with the unapodized mask. 

Given the apodized mask and the fiducial power spec¬ 
trum used to generate the simulations, we numerically 
evaluate the modified shape function using the 

expression in Eq. [18]. The 3} and 6j symbols are cal¬ 
culated using standard routines available in SLATEC nu¬ 
merical libray^. The modified shape functions evaluated 
and used in our analysis are depicted in Fig. [1]. 


C. Analysis and results 

In this section we test the recovery of Doppler boost 
vector from simulated data using the new estimator pre¬ 
sented in section IV B. We perform the analysis using the 
full sky simulations in addition to the analysis on masked 
maps to allow direct comparison. 


^ http;//healpix.sourceforge.net/ 

^ http;//irsa.ipac.caltech.edu/Missions/planck.html 
® http;//www.netlib.org/slatec/index.html 
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FIG. 2: The unapodized galactic mask used in the analysis is 
shown here. This {1,0} mask has an available sky fraction of 
fsky ~ 0.78, which remains practically same on apodization. 


For the masked analysis, we first evaluate the ensemble 
averaged BipoSH coefficients and their variance from the 
unboosted, noise added, masked simulations. While the 
average is used to subtract the biases due to the mask and 
spatially varying noise, the variance is used in the evalua¬ 
tion of the minimum variance estimator in Eq. [21]. The 
Doppler boost estimator is evaluated on each realization 
of Doppler boosted simulations, yielding an ensemble of 
estimates of the Doppler field P{n) = fi ■ n. While the 
direction of the Doppler boost is determined using the 
HEALPix subroutine remove_dipole, the Doppler boost 
amplitude is recovered by first estimating the power in 
the reconstructed Doppler field : |;5| = where 

/3l = Z^M l/^lMp/3. 

The Doppler power is not expected to vanish when es¬ 
timated from an unboosted sky, and this bias is termed 
reconstruction noise. In order to have an unbiased esti¬ 
mate of the Doppler amplitude the reconstruction noise 
needs to be subtracted from the estimated Doppler 
power. This bias can be dealt with semi-analytically in 
the full sky case (see Appendix A). For the masked anal¬ 
ysis we estimate the mean reconstruction noise by apply¬ 
ing the Doppler estimator on a set of masked unboosted 
simulations. The unbiased estimate of the Doppler boost 
amplitude is given by |/3| = 1.5y/(/3i — /3[^)/7r. The full 
sky analysis follows the same procedure as above, except 
that we don’t have to bias subtract the BipoSH coeffi¬ 
cients, since we only run the full sky analysis on simu¬ 
lations with isotropic noise. Note that one still needs to 
bias subtract the power in the Doppler field to arrive at 
the correct Doppler boost amplitude. 

Finally we run the analysis on two multipole ranges, 
[2,1000] and [2,2000]. At high multipoles I > 1000, the 
results can be potentially affected by point source masks. 
By demonstrating robustness against a range of multi¬ 
poles used in the analysis, we have shown that our results 
are not affected by the presence of point source masks. 

Results of the reconstruction of Doppler boost vector 
are summarized in Fig. [3]. Table II lists the mean re¬ 
covered Doppler amplitude and direction from different 
multipole bins. Note that the injected amplitude and 



p,/ii 


0.4 


0.3 


0.2 


0.1 


0 


Full sky/isotropic noise — 
Masked sky/isotropic noise — 
Masked sky/anisotropic noise — 




FIG. 3: This figure depicts the recovered Doppler boost pa¬ 
rameters from analysis on simulated GMB skies masked with 
the apodized Planck common mask. While the left column de¬ 
picts the histogram of the recovered Doppler power corrected 
for reconstruction noise bias, the right column depicts the his¬ 
togram of the cosine of angular separation /3o • /3, where /3o and 
/3 denotes the injected and recovered directions respectively. 
The top and bottom row denote the results from analysis in 
the multipole bins I = [2,1000] and [2, 2000] respectively. The 
injected Doppler power is denoted by a vertical grey line. 


Mask 

Al 

|/3| (xlO-D 

b° 

1 ° 

Full sky 

[2,1000] 

1.23 ± 0.33 

46 ± 14 

266 ± 26 

(isotropic noise) 

[2.2000] 

1.23 ± 0.18 

47.7 ± 8.2 

266 ± 13 

Common Mask 

[2,1000] 

1.28 ± 0.38 

45 ± 15 

262 ± 32 

(isotropic noise) 

[2.2000] 

1.28 ± 0.20 

46.6 ± 9.2 

262 ± 14 

Common mask 

[2,1000] 

1.27 ± 0.38 

44 ± 15 

263 ± 30 

(anisotropic noise) 

[2.2000] 

1.28 ± 0.21 

46.2 ± 9.2 

262 ± 15 


TABLE II: The mean recovered Doppler amplitude and direc¬ 
tion obtained from 1000 simulations are tabulated here. The 
boosted simulations were generated with a Doppler boost am¬ 
plitude |{?| = 1.23 X 10~^ pointing towards the galactic coordi¬ 
nates {I, b) — (264°, 48°). The errors quoted are the standard 
deviation computed from probability density function inferred 
from the ensemble of reconstructed Doppler power and direc¬ 
tion (l,b). 


direction are consistently recovered from the masked as 
well as full sky simulations with isotropic and anisotropic 
noise. We find that the error on the recovered Doppler 
amplitude and direction is larger when using smaller mul¬ 
tipole range I G [2,1000] as compared to estimates from 
using the larger multipole range I G [2, 2000] as expected. 
Similarly the error on estimates from masked sky is larger 
than the error on estimate from full sky, as expected due 
to the reduced sky fraction owing to the mask. Finally 
we reiterate that we don’t have to perform any fsky cor¬ 
rections to the Doppler amplitude as is usually required 
by other methods. 
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VI. CONCLUSIONS 

For any reliable cosmological analysis, masking the 
CMB maps is inevitable to avoid biases due to galac¬ 
tic foregrounds. The effects of masking on the angular 
power spectrum are easy to reverse using the MASTER 
algorithm [23]. 

Here we developed a formalism similar to the MAS¬ 
TER algorithm but extended to all BipoSH coefficients 
which are a generalization of the well known angular 
power spectrum. Using this formalism, we derived an 
expression which describes how masking modifies an ar¬ 
bitrary isotropy violating signal. Though we write down 
the mask coupling matrix for the BipoSH coefficients, we 
don’t try to invert the equation as done in the MASTER 
algorithm. We simplified the equation using the sym¬ 
metry properties of the mask and properties of Clebsch- 
Gordon coefficients. 

As a very specific example, we studied how the Doppler 
boost signal is modified due to the presence of the mask. 
We showed that by making reasonable assumptions it 
is possible to cast the BipoSH coefficients of a masked 
Doppler boosted sky in a form similar to the BipoSH 


coefficients of Doppler boosted full CMB sky, where the 
full sky shape function is replaced by the modified shape 
function. Eventually we obtained an estimator which 
recovers the Doppler signal from a masked CMB sky. 
This estimator naturally accounts for the loss of power 
due to masking, and hence does not require an additional 
fsky correction as required by other quadratic estimators. 

We generated a set of simulations whose statistical 
properties are tailored to match the 217 GHz Planck 
maps. We used an apodized form of Planck 2015 common 
mask for our analysis. Finally we evaluated the newly de¬ 
rived estimator on masked CMB skies to demonstrate it’s 
unbiased recovery of the injected signal. 

Finally we note that the method developed here are 
generic and can be used to reconstruct any arbitrary form 
of SI violation. 
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Appendix A: Full sky statistic 


Many isotropy violating signals can be written in the 
generic form [7], 


iLM _ aLM 
^hh - ^hh ' 




h ’ 


(Al) 


where denote the BipoSH coefficients measured 

from the data, denote the BipoSH coefficients of the 
statistically isotropic CMB sky, Plm denote the parame¬ 
ters of the isotropy violating field, and denotes the 
spectral shape function induced by the particular form 
of isotropy violation. 

For a statistically isotropic CMB field, = 0 for 

L > 0. Following Eq. [Al], we can define a naive estima¬ 
tor for the parameters /3lm as follows. 


1 JLM 


'^hh q/2 ^hh 


(A2) 


However, the above estimator is not optimized to re¬ 
duce the variance due to noise. It is possible to recover a 





signal with better signal to noise ratio by redefining the 
estimator as, 


aLM 

hh 'I'l 


(A3) 


where the weights be chosen so as to minimize 

the quantity, 


= 


{PlmPlm) ) 

, 2Cl> Cv aa 


= A^L + Cf 


(A4a) 

(A4b) 

(A4c) 


where fV/, is the reconstruction noise that is minimized by 
choosing appropriate weights, and denote the 

reconstructed and true power in the multipole L respec¬ 
tively, of any isotropy violating field /3(h). Minimizing 
the reconstruction noise results in down weighting the 
noisy modes while giving more weightage to the signal 
dominated modes. In arriving at Eq. [A4b] we used the 
covariance of BipoSH coefficients given by [4] , 

MM' 5 

(AS) 

which is valid when L ^ 0 and L + li + I 2 is even. 

This minimization of power is evaluated subject to the 
constraint resulting in a constrained min¬ 

imization problem, which is solved using the method of 
Lagrange multipliers. The weights are determined by 
minimizing the function. 


iKTAvlif) = 


L = C^/ -a 


.hh 


h' i> 
‘■1^2 


- 1 


(A6) 


where a is the Lagrange multiplier. Note that the weights 
for each BipoSH mode are assumed to be independent of 
each other in the following sense. 


dw 


hh 


dwh',, 




(A7) 


On setting the derivative of C with respect to the weights 
wp i, to zero and simplifying the resultant equation it can 
be shown that the weights that minimize the reconstruc¬ 
tion noise are given by, 


Whl 


CiM. 


E 

-hh 




(AS) 


Finally the reconstruction noise of the minimum variance 
estimator is given by, 


Nl 


ST' 

k 


(A9) 


Appendix B: Partial sky statistic 


We follow a procedure similar to the one described in 
Appendix A, to arrive at a minimum variance estimator 
in the presence of a mask. Following Eq. [19], we can 
define a minimum variance estimator for any isotropy 
violating field /3(h) as, 


/3lm = V 


LM 


ala 

■^hl 


hh 


hh t^LM 
^^LMh h 


(Bl) 


where .4^^ = — (A^^), A;^^ are the BipoSH co¬ 

efficients measured from the data, and (A^^) denotes 
the bias due to all known systematic effects in the data 
like masking, anisotropic noise, non-circular beam, etc., 
estimated from simulations which incorporate all these 
systematic effects. It is non-trivial to evaluate the co- 
variance of the BipoSH coefficients analytically, as 
can be done in the case of ideal full sky CMB. So we 
approximate the BipoSH covariance to be diagonal i.e.. 


{Af^j^A^jf )sim. « hhfh ^IiI[^121 '2^LL'5 mM' , (B2) 


where is the variance of BipoSH coefficients 

which is estimated from simulations. Note that the co- 
variance of BipoSH coefficients have an explicit M depen¬ 
dence which was absent in the full sky BipoSH covariance 
given in Eq. [A5]. We also assume that the weights for the 
different BipoSH modes are independent of each other as 
in Eq. [A7]. Using these approximations and following 
the same procedure discussed in Appendix A, one finds 
that the weights that minimize the variance of the esti¬ 
mator are given by. 


K 12 = 


E 


M 


E- 


^LM 

l' 

hh 

t^LM 
LMli li- 


(B3) 


Note that even though the weights themselves do not 
have an explicit M dependence, the effective weighting of 
the BipoSH coefficients inferred from data is M depen¬ 
dent owing to the presence of the modified shape function 
(see Eq. [Bl]). 





















